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The effects of dipolar interactions and molecular flexibility on the structure and phase behavior of 
bent-core molecular fluids are studied using Monte Carlo computer simulations. Some calculations 
of flexoelectric coefficients are also reported. The rigid cores of the model molecules consist of either 
five or seven soft spheres arranged in a 'V shape with external bend angle 7. With purely repulsive 
sphere-sphere interactions and 7 = 0° (linear molecules) the seven-sphere model exhibits isotropic, 
uniaxial nematic, smectic-A, and tilted phases. With 7 > 20° the smectic-A phase disappears, while 
the system with 7 > 40° shows a direct tilted smectic-isotropic fluid transition. The addition of 
electrostatic interactions between transverse dipole moments on the apical spheres is generally seen 
to reduce the degree of tilt in the smectic and solid phases, destabilize the nematic and smectic- 
A phases of linear molecules, and destabilize the tilted smectic-B phase of bent-core molecules. 
The effects of adding three-segment flexible tails to the ends of five-sphere bent-core molecules are 
examined using configurational-bias Monte Carlo simulations. Only isotropic and smectic phases 
are observed. On the one hand, molecular flexibility gives rise to pronounced fluctuations in the 
smectic-layer structure, bringing the simulated system in better correspondence with real materials; 
on the other hand, the smectic phase shows almost no tilt. Lastly, the flexoelectric coefficients of 
various nematic phases - with and without attractive sphere-sphere interactions - are presented. 
The results are encouraging, but the computational effort required is a drawback associated with 
the use of fluctuation relations. 



I. INTRODUCTION 

There has been a surging interest in bent-core liquid 
crystals since their discovery in 1996 k Typically, these 
materials consist of molecules comprising rigid, banana- 
shaped cores made up of a conjugated system of linked 
aryl groups, and flexible alkyl or alkoxy tails attached to 
each end. The molecules are usually achiral and possess 
electric dipole moments parallel with the molecular C2 
axes. One of the most intriguing properties of these com- 
pounds is that, in some cases, chiral ferroelectric or anti- 
ferroelectric smectic phases can be observedi*^. The chi- 
rality arises because the molecules tilt within the smec- 
tic layers; in chiral (anti)ferroelectric phases all of the 
molecules tilt in the same sense with respect to the layer 
polarization vector. It is not yet clear what is responsi- 
ble for this spontaneous symmetry-breaking process, al- 
though a variety of explanations has been proposed. One 
popular explanation involves the long-range dipole-dipolc 
interactional, while recent theoretical work has identified 
a central role for dispersion interactions 5 . Other possible 
explanations include entropic 'free-volume' mechanisms 6 
in which an antiferroelectric ordering of the smectic-layer 
polarizations affords more room for layer fluctuations, 
and mechanisms in which the molecules themselves spon- 
taneously select chiral molecular conformations^. 

There is a growing simulation literature on bent- 
core liquid crystals. One of the most simple bent- 
core molecular models is a dimer made up of two hard 
spherocylinders^A This system exhibits isotropic, ne- 
matic, smectic, and crystalline phases, but no tilted 
phases. The Gay-Berne dimer model has been studied 



extensively, with and without molecular dipoles. In the 
works by Memmer- and Johnston et o/^S*ii, isotropic, 
nematic, tilted smectic, and helical phases were found, 
depending on the molecular bend angleS*^ and the mag- 
nitude of the dipole momenljAi. Xu et al. studied com- 
posite molecules made up of repulsive soft spheres, and 
found isotropic and tilted crystalline phases^. More re- 
cently, we have studied composite molecules made up of 
Lennard-Jones spheres - so called 'composite Lennard- 
Jones molecules' (CLJMs) - which exhibit isotropic, ne- 
matic, tilted smectic, and tilted crystalline phases^. 

In this work the effects of molecular dipole moments 
and molecular flexibility on the phase behavior of model 
bent-core molecules are studied using computer simula- 
tions. The model (to be detailed in Section [H} con- 
sists of a rigid 'V'-shaped core of soft spheres with a 
point dipole moment oriented along the C2 axis. Molec- 
ular flexibility is included by the addition of short flexi- 
ble tails to either end. There is a substantial literature 
on the effects of these molecular characteristics on lin- 
ear molecules. In hard-spherocylinder fluids, the addi- 
tion of longitudinal molecular dipoles is seen to destabi- 
lize the nematic phase, and can even destabilize smectic 
phases if the dipoles are displaced toward the ends of the 
moleculesiii 5 -; transverse dipoles also destabilize the ne- 
matic phase with respect to the smectic Ai&. Gay-Berne 
ellipsoids with longitudinal point dipoles show a stabiliza- 
tion of the nematic phase with respect to the isotropic 
phase as the dipoles are moved from the centers of the 
molecules to the endsii, and can exhibit antiferroelectric 
smectic phases with striped structures^ 8 . Tilted polar 
smectic phases have been reported in fluids of Gay-Berne 
molecules with transverse dipole momentaiS,. As far as 



molecular flexibility is concerned, the general consensus 
is that the introduction of flexible tail groups destabi- 
lizes the nematic phases of hard spherocylinders 2 ^, fused 
hard-sphere chains2ii22, Gay-Berne^ 3 ., and soft-sphere 
chains 24 . Interestingly, the simultaneous presence of flex- 
ible tails and molecular dipole moments can lead to a sta- 
bilization of the nematic phas o 2 ^ 2 ^ . With regard to bent- 
core molecules, Johnston et al. have shown that the pres- 
ence of transverse molecular dipoles on Gay-Berne dimers 
stabilizes the smectic phases at the expense of nematic 
phases, increases the tilt angle in tilted smectic phases, 
and can induce long-range polar ordering^. In the cur- 
rent work we will show that for the bent-core soft-sphere 
models considered, the additions of dipolar interactions 
and flexible tails both destabilize the nematic phase, and 
that the dipolar interactions reduce the degree of molec- 
ular tilt in smectic phases. 
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FIG. 1: The molecular models studied in this work: (a) the 
CSSM model; (b) the CSSMT model. Also shown are the 
molecular axes, a, b, and c. 



II. MOLECULAR MODEL AND SIMULATION 
METHODS 



This paper also reports our attempts to measure 
the flexoelectric coefficients^^ of model bent-core 
molecules. There are relatively few accounts of such 
measurements in the literature. Experimentally, the de- 
termination of these quantities is highly non-trivialS2^, 
mainly due to the fact that the flexoelectric coefficients 
are not measured directly, but rather in linear combina- 
tions or as ratios involving elastic constants. In simula- 
tions, the flexoelectric coefficients of pear-shaped Gay- 
Berne ellipsoid/Lennard- Jones sphere molecules have 
been measured directly using expressions involving the 
direct correlation functional. The coefficients for a 
similar model were studied in simulations using fluc- 
tuation expressions^. The compound 5CB has been 
studied using a parameterized dipolar Gay-Berne model 
and the Percus-Yevick closure of the Ornstein-Zcrnike 
equationS^ii. The flexoelectric coefficients were com- 
puted using the direct correlation function route, and 
the results compared moderately well with experiments^. 
A very recent simulation study of PCH5 - using a fully 
atomistic molecular model - employed fluctuation for- 
mulae which yielded results in good agreement with 
experiment 3 -. In the present work we critically assess 
the reliability of the fluctuation route in the context of 
our model systems, and show that the bend flexoelectric 
coefficients for non-polar molecules can be comparable to 
those measured in experiments, reflecting the significant 
role of molecular packing in dense liquids. 



This paper is organized as follows. In Section [H] the 
molecular model to be studied is fully defined, and the 
required simulation methods are described. Simulation 
results for rigid linear molecules are presented in Section 
IIII Al and those for rigid bent-core molecules in Sections 
IIII Bl and IIII CI The effects of molecular flexibility are 
considered in Section IIII Dl and flexoelectricity is dis- 
cussed in Section Till El Section llVl concludes the paper. 



The molecular models considered in this work are 
shown schematically in Fig. ^ The most basic model 
[Fig. IHa)] consists of a rigid core of seven soft spheres 
arranged in a 'V shape with external bond angle 7 de- 
fined such that ei ■ e2 = cos (180° —7), where ei and 
e 2 are unit vectors pointing along the two 'arms' of 
the molecule. The sphere-sphere interaction potential 
is taken to be the repulsive part of the Lennard- Jones 
(12,6), i.e., 

u ss (r)^4e(^y\ (1) 

where r is the sphere-sphere separation, e is an energy 
parameter, and a is the sphere 'diameter'. In this work 
the intramolecular sphere-sphere bond length is set equal 
to la. To identify the molecular axes that will be aligned 
in orientationally ordered phases, we define three unit 
vectors associated with the rigid cores of the molecules, 
as illustrated in Fig. ^ These vectors are given by a = 
(ei-e 2 )/|ei — e 2 |, b = (ei +e 2 )/|ei +e 2 | and c = aAb. 

Many real bent-core molecules possess a transverse 
electric dipole moment aligned along the C 2 molecular 
symmetry axis. To represent dipole-dipole interactions 
we have considered the model represented in Fig. ^a), 
where a point dipole moment is placed at the center of 
the apical sphere along the C 2 symmetry axis (b). The 
dipole-dipole interaction is, 

H1H2 3(/xi • r)(/i 2 • r) 
UddfaHuVh) = ^ L , (2) 

where r is the pair separation vector, r = |r|, and fit = 
/ibi is the dipole vector on particle i. For brevity, we 
shall refer to these molecules as 'composite soft-sphere 
molecules' (CSSMs). 

Common additional features of real bent-core 
molecules include alkyl or alkoxy chains (typically 
3-12 carbons in length) attached to both ends of each 
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molecule. To represent these tails we consider the 
addition of four extra spheres at each end of the model 
molecules, as represented in Fig. db). The resulting 
tail segments are allowed to pivot under the following 
constraints: the tail bond length is equal to 0.6c, which 
corresponds to the ratio of the carbon-carbon bond 
length to the diameter of the aromatic ring in real 
bent-core liquid crystals; each tail segment is oriented 
at the tetrahedral angle, cos -1 (—5) = 109.47°, with 
respect to its neighboring segments, mimicking the 
bond angles in a simple hydrocarbon tail. Clearly 
this extension will cause a considerable increase in the 
number of interactions to be evaluated, so to make 
the simulations tractable we have made the following 
simplifications. Firstly, we have reduced the size of 
the rigid bent-core to five spheres. (With the addition 
of the tail segments, this actually makes the effective 
elongations for all of the models considered more com- 
parable.) Secondly, we neglect interactions between tails 
on different molecules. Within the Lorentz-Bcrthclot 
mixing rules^, this corresponds to setting the tail-sphere 
diameter to zero. Denoting the core and tail spheres by 
'c' and 't', respectively, the interaction potential is, 

<i(r)=4e(^) 12 , (3) 

with i,j — c or t, a cc = a, a u = 0, a ct = \(o cc + 
c«) = \<J, and for simplicity the energy parameter e 
is the same for all pairs. This is clearly a very crude 
representation of molecular flexibility, but it has proven 
to be an appropriate means of extending the range of 
applicability of simple liquid-crystal models^fi^. Even 
though this is a simple model, configurational-bias MC 
techniques are required to simulate the system efficiently; 
these are summarized in Section III Al For brevity, we 
shall refer to these molecules as 'composite soft-sphere 
molecules with tails' (CSSMTs). 

Reduced units for these systems arc defined as follows: 
reduced molecular density, p* = N m a 3 /V, where N m is 
the number of molecules and V is the system volume; 
reduced temperature, T* = ksT/e; reduced pressure, 
p* = p<r 3 /e; reduced dipole moment, p* = \J fi 2 /ea 3 . 

A. Monte Carlo 

The phase behavior of the model systems was inves- 
tigated using constant-pressure (NPT) and constant- 
volume (NVT) Metropolis MC simulations^. In all 
of the simulations reported in this work, the number 
of molecules was N m = 400, with initial high-density 
crystalline configurations consisting of four layers of 100 
molecules. The general approach was to equilibrate the 
system at low temperature (T* ~ 1) and high density 
(p* ~ 0.14) using NVT simulations in a cuboidal simu- 
lation cell with dimensions L x — L y L z and volume 
V = L x L y L z , and then switching over to NPT simula- 
tions at a fixed pressure of p* = 4; earlier work indicated 



that all of the expected liquid-crystalline phases could be 
stabilized in a CLJM system at this pressure^. Simula- 
tions along the isobar were carried out at progressively 
higher temperatures in order to locate phase transitions 
between solid, smectic, nematic and isotropic phases; 
transitions were identified by measuring the equation of 
state (density as a function of temperature) and relevant 
order parameters (detailed below). Cooling runs were 
carried out using cuboidal and/or cubic simulation cells 
to confirm the existence and nature of the transitions. 
As explained in Ref. Il.'jl spot checks on the stress tensor 
showed that the cuboidal/cubic cells did not mechani- 
cally destabilize liquid-crystalline phases. 

For the CSSM system, straightforward simulation tech- 
niques were employed as detailed in earlier worki^; single- 
particle translation and rotation moves, and volume 
moves (in In V) , were generated with respective maxi- 
mum displacement parameters to achieve ~ 50% accep- 
tance rates. The long-range dipolar interactions were 
handled using Ewald summations with conducting ('tin- 
foil') boundary conditions^. 

For the CSSMT system the CBMC techni que wa s im- 
plemented and optimized as described in Refs. 37 38. Tail 
conformations were sampled by generating 5 trial orien- 
tations per segment per MC move. Translational and ro- 
tational displacement parameters were adjusted to give 
an acceptance ratio of 10%; in NPT simulations the vol- 
ume moves were adjusted to give a 50% acceptance ratio. 
The computational effort required to simulate this system 
was considerable. To carry out a MC sweep consisting 
of one attempted translation and rotation per molecule, 
and one volume move, took approximately 4 seconds on 
a 2.2 GHz Intel Xeon processor; to achieve equilibration 
at each state point required at least 10 5 MC sweeps 

To monitor orientational order, the order tensors 

Qaa = \ ^2iLi(3 aa ~~ 1) f° r eacri 01 the molecular axes 
a = a, b, and c were diagonalized yielding the eigenval- 
ues A~ < A° < A+, and the corresponding orthonormal 
eigenvectors, n~, n° , and n+— . The molecular x, y, and 
z axes were then assigned in order of increasing A„, and 
the laboratory axes, X, Y, and Z, were identified with 
the corresponding directors. In practice this almost in- 
variably meant that the molecular z axis was the 'long' 
axis a and the x and y axes were b and c, and that the 
laboratory z axis was n+ . The usual nematic and biaxial 
order parameters - S and Q\ 2 - are then given by, 

S = Z Q zz Z, (4) 

Q22 = Jj(^ ' Qza; ' X + Y • Qyy • Y 

— X • Qyy • X — Y • Q xx ■ Y). (5) 

In a perfect uniaxial nematic phase, S = 1 and Q22 = 0> 
whereas in a perfect biaxial phase, S = 1 and Q22 = !• 
In practice no biaxial ordering was detected in any of 
the simulations, and so we will not report the numerical 
values of Q\ 2 (which are less than ~ 0.1). Some ad- 
ditional measured observables include the polarization, 
P = fj, X^bj, an d the intermolecular torque-density ten- 
sor, II. In particular, these quantities are required for 
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the calculation of flexoelectric coefficients, full details of 
which will be given in Section llll El 



III. RESULTS 
A. CSSMs with 7 = 0° 

The equation of state and order parameters for apo- 
lar (/i* = 0) linear CSSMs along an isobar with p* = 4 
are shown in Figs.EJa) andEJb), respectively. In order 
of increasing temperature we find a high-density tilted 
phase (0.5 < T* < 1.5), a smcctic A (2.0 < T* < 2.5), 
a uniaxial nematic (2.75 < T* < 4), and ultimately the 
isotropic phase (T* > 4.5). Simulation snapshots are 
shown in Fig. [3] In the tilted phase, the molecules are 
arranged in layers tilted by about 60° with respect to 
each other, as shown in Fig. a). This 'herringbone' 
structure clearly allows close-packing of the constituent 
spheres. It is difficult to resolve the molecules in to layers 
unambiguously, but it is quite clear from examining sim- 
ulation snapshots that there is no long-range crystalline 
order. In the absence of such order we therefore classify 
this phase as tilted smectic B, if only to indicate that it 
is not crystalline 28 . We found no stable crystalline phase 
at T* > 0.5. Referring to Fig.^b), the jump in nematic 
order parameter in the temperature range 1.5 < T* < 2.0 
is due to the transition from the tilted smectic-B phase 
to the untilted smectic-A phase; the drop in the range 
4 < T* < 4.5 signals the smectic A-nematic transition. 
In all, the results for this system are in good qualita- 
tive correspondence with those for a whole host of similar 
(linear) molecular models, including soft-sphere chains^, 
and Lennard-Jones chains^i. The results are also compa- 
rable to those presented in Ref.^^for seven-sphere CLJM 
fluids in which the sphere-sphere interaction is given by 
4e[(cr/r) 12 - (cr/r) 6 ]. Qualitatively, the CSSM and CLJM 
systems are very similar, but in the latter case the phase 
transitions are shifted to higher temperatures due to the 
attractive component of the interaction potential. 

With the addition of a small molecular dipole (/1* = 1) 
we see little qualitative difference in the equation of state 
at p* = 4, as shown in Fig.^c). Despite the small change 
in the equation of state, the gross structure of the tilted 
smectic-B phase is quite different from that in the apolar 
system. The low-temperature smectic-B phase is not so 
strongly tilted as in the apolar system, exhibiting a tilt 
angle with respect to the layer normal of ~ 20°. This is 
most likely to allow dipoles on neighboring molecules to 
attain the low-energy 'nose-to-tail' conformation within 
the plane of the layer. The nematic order parameter - 
shown in Fig. - is relatively high at temperatures 
T* < 1.5 due to the reduced degree of tilt. The tilted 
smectic B-smectic A and smectic A-nematic transitions 
are signaled by changes in S at 1.5 < T* < 2 and 4 < 
T* < 4.5, respectively. 

With a large dipole moment (p* — 2) and p* = 4 we 
see some dramatic differences in the phase behavior: the 
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FIG. 2: Equations of state [(a),(c),(e)] and order parameters 
[(b),(d),(f)] for CSSM systems with 7 = 0° along an isobar 
with p* = 4: (a),(b) p* = 0; (c),(d) p* = 1; (e)m(f) p* = 2. 
The symbols denote different phases: solid (down triangles); 
tilted smectic B (up triangles); smectic A (diamonds); ne- 
matic (squares); isotropic (circles). 





FIG. 3: Simulation snapshots of the CSSM system with 7 = 
0° (linear molecules) along an isobar with p* = 4: (a) tilted 
smectic-B (p* = 0, T* = 2); (b) smectic-A (p* = 0, T* = 
2.5); (c) nematic (p* — 0, T* =4); (d) tilted smectic-B 
(A** = 1. T* = 2). 



nematic and smectic-A phases are completely absent, and 
there is instead a distinct transition between two high- 
density layered phases in the range 2 < T* < 2.5. Simu- 
lation snapshots at temperatures of T* = 2 and T* = 4 
are shown in Fig. Examination of the layers [Fig.^Jc) 
and^d)] shows that at T* < 2 the system is in a crys- 
talline phase, with apparently long-range positional order 
within the layers. At T* = 4 the in-layer ordering is qual- 
itatively different, showing short-range positional corre- 
lations and defects that destroy long-range positional or- 
der. The equation of state and nematic order parameter 
are shown in Figs.HJe) and|2If). We assign the branches 
in the equation of state as corresponding to crystalline 
(0.5 < T* < 2.0), tilted smectic-B (2.5 <T*< 6.0), and 
isotropic (T* > 6.5) phases. 

To summarize, the addition of dipolar interactions to 
linear seven-sphere molecules leads to a reduction in the 
degree of tilt in the low-temperature smectic-B phase. 
With high dipole moments, the smectic-A and nematic 
phases disappear, and a solid-smectic B transition is 
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FIG. 4: Simulation snapshots of the CSSM system with 7 = 
0° and /1* = 2: (a) crystalline phase (T* = 2); (b) tilted 
smectic-B phase (T* = 4) ; (c) a typical layer in the crystalline 
phase (T* = 2); a typical layer in the tilted smectic-B phase 
(T* =4). In (c) and (d) the short black lines indicate the 
orientations of the molecular dipole moments. 
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FIG. 5: Equations of state [(a),(c),(e)] and order parameters 
[(b),(d),(f)] for CSSM systems with 7 = 20° and 7 = 40° 
along an isobar with p* = 4: (a),(b) 7 = 20° and fi* = 0; 
(c),(d) 7 = 20° and fi* = 1; (e),(f) 7 = 40° and fi* = 0. 
The symbols denote different phases: tilted smectic B (up 
triangles); nematic (squares); isotropic (circles). 



shifted in to the range of temperatures considered in this 
work. 




B. CSSMs with 7 = 20° 

The equations of state and nematic order parameters 
for CSSM systems with 7 = 20° along an isobar with 
p* = 4 are shown in Fig. El Results are shown for two 
dipole moments, fi* = [Figs. Eta) and Etb)] and fj,* = 1 
[Figs. Etc) and Eld)]. Both the apolar and polar systems 
exhibit tilted smectic-B, nematic, and isotropic phases; 
examples of the smectic and nematic phases in the fj,* = 
system are illustrated in Fig. El In the smectic-B phases 
it was observed that the degree of molecular tilt with re- 
spect to the layer normal is far greater in the apolar case 
(~ 53°) than in the polar case (< 20°). The smectic B- 
nematic phase transition appears to be more pronounced 
in the apolar system than in the polar system, as ev- 
idenced by the associated features in the equations of 
state and in the variations of the nematic order parame- 
ters. Some simulations were attempted with fi* = 2 but 
these suffered from convergence problems; simulations 
with different initial configurations failed to converge on 
to the same branch of the equation of state. It is pos- 
sible that this was due to the combination of the steric 
dipole (molecular bend) and the 'electric' dipole resulting 
in strong anisotropic interactions and prohibitively slow 
convergence. 

A comparison of Figs. [21 and El shows that the pres- 
ence of a modest molecular bend leads to the smectic-A 
phase being destabilized. This same trend was observed 
in simulations of the CLJM system^. The introduction 
of dipolar interactions to the bent-core model then seems 
to stabilize the nematic phase slightly in favor of the 
smectic B. The smectic-B phases themselves are tilted, 
but the degree of tilt is reduced significantly upon the 
addition of dipolar interactions. This perhaps provides 




FIG. 6: Simulation snapshots of the CSSM system with 7 = 
20° along an isobar with p* = 4: (a) two views of the tilted 
smectic B phase at T* = 1; (b) nematic phase at T* = 3. 



a clue as to why dipolar interactions apparently disfavor 
the smectic B; the bent-cores want to form a tilted phase, 
but the dipolar interactions want an untilted phase as ex- 
plained in Section lill Al 



C. CSSMs with 7 = 40° 

The convergence problems encountered with the 7 = 
20° system were exacerbated by an increase of the bend 
angle to 7 = 40°. In this case it was only possible to 
achieve reliable results for the apolar system. Results 
at \i* = and p* = 4 are shown in Fig. EJe) and Elf). 
The equation of state and order parameters show only 
two branches, which correspond to tilted smectic-B and 
isotropic phases. This is very similar to the situation 
in the CLJM system with the same bend angled, albeit 
with the CSSM system undergoing a phase transition at 
lower temperature. 



D. CSSMTs 

We performed NPT simulations of CSSMT systems 
with bend angles of 7 = 0°, 20°, and 40° along an iso- 
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FIG. 7: Simulation snapshots of CSSMT systems at T* = 1 
and p* = 4: (a) 7 = 0°; (b) 7 = 20°; (c) 7 = 40°. 



bar with p* = 4. Due to the computational effort re- 
quired for these simulations we were not able to map 
out equations of state as comprehensive as those for the 
CSSM systems. The simulation results are presented in 
Table [I] The only phases observed in our simulations 
were smectic and isotropic. In all cases the smectics were 
stable at T* < 1.5, and the smectic-isotropic transition 
occurred in the range 1.5 < T* < 2. Some simulation 
snapshots of the smectic phases are shown in Fig. [7| In 
all cases the translational ordering within the layers was 
of the smectic-B type, i.e., local hexagonal coordination. 
The systems with 7 = 0° and 7 = 20° showed no ap- 
preciable molecular tilt within the smectic layers, while 
the 7 = 40° system showed an unusual 'grain-boundary' 
structure between clearly demarcated domains of untilted 
smectic. Larger-scale simulations will be required to de- 
termine whether this is a signal of a long-wavelength 
modulated structure. It is particularly striking that the 
smectic-layer fluctuations are much larger than those in 
the CSSM (and CLJM 13 ) systems. Indeed, one criticism 
of the latter models - and other rigid-rod models - is that 
the smectics are too well ordered. Unsurprisingly, the in- 
troduction of molecular flexibility has improved the cor- 
respondence between simulated smectic structures, and 
those inferred from light-scattering experiments on com- 
mon (flexible or semiflexible) mesogens 2 ^. 

The conformations of the flexible tail groups were in- 
vestigated using some simple measures. The extension 
of each tail was identified with the distance, Z, between 
the first and fourth joint (the black spheres in Fig. QJ. 
The probability density function, p(l), is shown in Fig. [S] 
for all of the CSSMT systems at T* = 1 (smectic B) and 
T* = 3 (isotropic) . Each function shows peaks at l/a = 1 
and l/a ~ 1.5. For a perfect cis conformation the tail ex- 
tension is (5/3) times the bond length, while for the trans 
conformation it is yl9/3 times the bond length. With 
the bond length being 0.6a, these distances correspond to 
l/a = 1 and l/a ~ 1.51, respectively. Firstly, the cis con- 
formation is clearly the more favorable, presumably be- 
cause the molecules strive to attain the shortest effective 
elongation to minimize excluded- volume interactions. In- 
terestingly, for each system the cis conformation appears 
slightly more favorable in the isotropic phase than in the 
smectic phase. This may be due to the opportunity for 
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FIG. 8: Tail-length distribution functions for the terminal 
four-center tails in CSSMT systems along an isobar with p* = 
1.0: (bottom) 7 = 0°; (middle) 7 = 20°; (top) 7 = 40°. 
The curves for each bond angle are displaced by 10 units 
along the ordinate for clarity. In each case the solid lines 
correspond to T* = 1 and the dashed lines to T* = 3. With 
bond lengths of 0.6<r the pure cis conformation corresponds 
to / = (5/3) x 0.6cr = a, and the pure trans conformation to 
I = y/19/3 x 0.6a ~ 1.51cr. 

interdigitation of the tails with the cores in neighboring 
smectic layers, which would explain the accompanying 
increase in the occurrence of the trans conformation. 

The effective molecular bend angle was determined by 
calculating the dot product of the two unit vectors linking 
the apical sphere to the terminal tail units, Ui and U2; 
the required angle is then 7 ff = cos -1 (— m • 112). Values 
of 7 c ff are reported in Table [I] For each system 7 e g > 7 
which shows that the tails must curl up in such a way 
as to make the molecule more banana shaped. In the 
smectic phases this can be observed directly in Fig. \7\ 
where the tails prefer to be oriented in the plane of the 
smectic layers, rather than pointing straight down toward 
the neighboring layers. This idea is confirmed by the fact 
that the effective molecular bend is more pronounced in 
the smectic phase than in the isotropic phase. 

In summary, the addition of molecular flexibility re- 
sults in the disappearance of the nematic phase, and in 
the case of the linear molecules, the smectic-A phase 
as well. This is in good qualitative agreement with 
the trends observed in a variety of other liquid-crystal 
modela 20 i 21 i 22 i 23 i 24 i 26 . The smectic-layer structures of the 
flexible model systems correspond more closely to those 
in real smectic liquid crystals. 



E. Flexoelectric coefficients 

In 1969 Meyer predicted the existence of what is now 
known as the flexoelectric effect in nematic liquid crys- 
tals, in which long-wavelength distortions of the local 
molecular alignment (director field, n) give rise to a bulk 
polarization, P 27 . The textbook explanation of the effect 
is that if the director field possesses curvature then an 
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TABLE I: Results from NPT simulations of the CSSMT system along an isobar with p* — 4. Digits in brackets denote the 
estimated statistical uncertainty in the last figure. 



7 / degrees 


T* 


P* 


S 


I /a 


7cff / degrees 





1.0 


0.1098(4) 


0.918 


5.48(1) 


48.0(3) 





1.5 


0.1018(4) 


0.725 


5.82(1) 


29.6(2) 





2.0 


0.1067(7) 


0.091 


5.40(3) 


31.8(3) 





2.5 


0.1019(7) 


0.097 


5.45(3) 


32.3(3) 





3.0 


0.0981(7) 


0.063 


5.44(2) 


32.5(3) 


20 


1.0 


0.1128(2) 


0.865 


5.41(1) 


50.9(4) 


20 


1.5 


0.1101(2) 


0.787 


5.51(1) 


43.3(2) 


20 


2.0 


0.1057(7) 


0.200 


5.35(3) 


37.0(5) 


20 


2.5 


0.1010(3) 


0.062 


5.34(3) 


37.4(4) 


20 


.3.0 


0.0976(7) 


0.054 


5.34(2) 


37.6(4) 


40 


1.0 


0.1161(3) 


0.879 


4.98(1) 


65.2(4) 


40 


1.5 


0.1084(1) 


0.833 


4.93(1) 


60.7(2) 


40 


2.0 


0.1048(6) 


0.067 


5.08(3) 


49.5(3) 


40 


2.5 


0.1008(6) 


0.081 


5.10(2) 


49.2(6) 


40 


3.0 


0.0974(7) 


0.057 


5.11(2) 


49.4(5) 



asymmetric molecular shape can dictate a favorable lo- 
cal packing arrangement which, in the presence of molec- 
ular dipoles, may give rise to a polarization^. With a 
splay deformation (V • n ^ 0) wedge-shaped molecules 
with longitudinal dipole moments pack most efficiently 
when the dipoles are aligned. With a bend deformation 
(|V A n 7^ 0) banana-shaped molecules with transverse 
dipole moments are arranged preferentially to give a net 
polarization. General symmetry arguments lead to the 
following relationship between the polarization density, 
p = V~ 1 P (units C m~ 2 ), and the lowest order deforma- 
tions of the director field, 

p = ei(V-n)n + e 3 (VAn) An. (6) 

where e\ and are the splay and bend flexoelectric 
coefficients, respectively, with units of C m _1 . Histor- 
ically there is some ambiguity in the sign of to be 
clear, throughout this work we employ the convention 
used by Meyer in his original studySi, Nemtsov and Os- 
ipov in their analysis of flexoelectricity in the context 
of linear-response theory^, and de Gennes and Prost in 
their canonical texljSS. Allen and Masters have supplied 
a comprehensive account of various simulation methods 
for measuring the flexoelectric coefficients^. Following 
the sign conventions in Ref.^3 we have calculated e\ and 
e-i using the relationships, 

6! = \$V~ X {{P z U xy ) - (P Z H yx )) , (7) 

e 3 = \PV~ X ({P y U zx ) - (P x Il zy )) . (8) 

where II = — '^2 i< j ^ijfij is the orientational stress den- 
sity tensor, ry = — Vj is the intermolecular separa- 
tion vector, and Ty is the torque on molecule i due to 
molecule j. Tij was calculated as a sum of moments of 
the sphere-sphere interactions about the apical sphere, 
and all vectors and tensors were calculated in a frame in 
which the laboratory z axis coincides with the nematic 



director, n+. ft is easy to show that the combinations 
(P z H xy — P z H- yx ) and (P y Tl zx — P x H zy ) are invariant with 
respect to a rotation of the x and y axes about n, and 
so the assignments of the x and y axes are arbitrary, ft 
should be noted that we have only calculated the flex- 
oelectric coefficients for non-polar systems. The 'steric 
dipole' is still parallel to b in Fig. ^ but there are no 
electrostatic dipole-dipole interactions. The polarization 
is given by P = /1 ^ii an d carries the trivial factor 
of \x by virtue of there being no electrostatic interactions. 

The flexoelectric coefficients have been calculated in 
the nematic phases of non-polar CSSM and CLJMi£ sys- 
tems as a function of the molecular bend angle, 7. For 
the purposes of comparison, the CSSM system has been 
studied at a fixed density and temperature for which 
the nematic phase is stable at several values of 7. An 
examination of Figs. EI a) and EJa) shows that in the 
range 0° < 7 < 20°, the nematic phase is stable at 
temperatures and densities in the regions of T* ^3 
and p* ~ 0.085, respectively. A particular state point 
from the 7 = 0° system was selected arbitrarily for all 
of the simulations, this being T* = 3 and p* — 0.0888. 
Canonical (NVT) simulations were used to equilibrate 
nematic phases for systems with bend angles in the range 
0° < 7 < 25°. The nematic order parameter, S, is shown 
as a function of 7 in Fig. G2a). We found that S could 
be fitted with a power law, 

S( 7 )=S(0)(l-£) a , (9) 

where S'(O) is the order parameter for linear molecules, 
7 C is a critical bend angle above which the nematic phase 
is no longer thermodynamically stable, and a is a specific 
exponent. The fit is shown in Fig.[!Ja); the fit parameters 
were 5(0) = 0.881(4), 7c = 25.6(2)°, and a = 0.104(7). 
A similar procedure was carried out for the CLJM system 
studied in Ref. 13J. Nematic phases were simulated at 



8 







■ M 




■ (o) 










: ^'^^c : 


■ (d) W 



5 10 15 20 25 30 5 10 15 20 25 30 35 

y/degrees y/degrees 



FIG. 9: Nematic order parameters [(a) and (c)] and flex- 
oelectric coefficients [(b) and (d)] for the CSSM system at 
p* = 0.0888 and T* = 3 [(a) and (b)] and the CLJM system 
at p* = 0.1155 and T* = 6.5 [(c) and (d)]. In (b) and (d) 
the flexoelectric coefficients are those corresponding to splay 
deformations (filled symbols) and bend deformations (open 
symbols). The lines are spline fits to guide the eye. 



T* = 6.5 and p* = 0.1155 for systems with 0° < 7 < 30°. 
The nematic order parameter is shown as a function of 
7 in Fig. Ifflc). A power-law fit yielded the parameters 
5(0) = 0.961(4), 7c = 29.6(1)°, and a = 0.143(4); the fit 
is shown in Fig. |Sfc). 

Equations Q and (JSJ) were evaluated using results 
from immense NVT-MC simulations. In general we car- 
ried out runs consisting of ~ 4 x 10 6 attempted MC 
translations and rotations per molecule. In all cases the 
components of P fluctuated about zero, but long runs 
were required to ensure that the average polarization, 
(P), was almost zero. Results for the reduced flexo- 
electric coefficients = ei<r 2 /fi and = e 3 a 2 /p in 
the CSSM and CLJM systems are presented in Table 
Hll and Fig. Also included in Table [D] are the com- 
ponents of the average polarization, P, and estimated 
uncertainties. In both the CSSM and CLJM systems, 
the measured flexoelectric coefficients are small when the 
molecular bend angles are less than about 20°. For more 
pronounced bend angles, both the splay and the bend 
coefficients deviate significantly from zero. The results 
for 63 are encouraging, but there is some concern over 
the measured values of e\; naively we would expect that 
the splay coefficients should be small for banana-shaped 
molecules with transverse dipoles. Table [H] yields a valu- 
able clue; the measured flexoelectric coefficients deviate 
significantly from zero in those simulations where there 
is a more pronounced average polarization. Therefore, it 
may be that the simulations are still not long enough to 
ensure a reliable evaluation of the fluctuation formulae. 
It is perhaps worth pointing out that during the course 
of the atomistic simulations performed by Cheung et al. 
there was a 'small net polarization—: the magnitudes of 
the average polarizations (P ~ 10 -28 C m^ 1 ) and the 
average molecular dipole moments (/x ~ 10~ 29 C m _1 ) 
give reduced polarizations P/fi ~ 10, which are large 
compared to the average polarizations reported in Table 



Hll It is possible that neither the estimates of Cheung 
et al. nor the current estimates are particularly reliable. 
We attempted to make ad hoc corrections to the fluctu- 
ation formulae in Eqns. J7J and JSJl by evaluating terms 
like ((P a — (P Q ))(II / 3 7 — (Il^y))), but these resulted in 
insignificant changes to the values of ei and e^. 

Notwithstanding the potential problems highlighted 
above, we can attempt to make some useful comments 
on the measured values of 63. At the highest bend an- 
gles, 7 = 20°-30°, the magnitude of the reduced bend 
flexoelectric coefficient is in the region of 0.5. With typ- 
ical values of a ~ 0.5 nm and p ~ 1 Dcbye, this reduced 
value corresponds to a real bend flexoelectric coefficient 
of e3 ~ 7 pC m _1 . This is in good agreement with typical 
values of e3 ~ 10 pC m _1 measured in experiments 2 ^!!. 
It is therefore reasonable to suggest that the steric or 
packing contributions to the flexoelectric effect are sig- 
nificant. The roles of dipolar and quadrupolar electro- 
static interactions must surely be at least as significant, 
and these should be studied in a systematic fashion. For 
now, though, we conclude that short-range interactions 
between bent-core molecules are as important in giving 
rise to flexoelectricity as they are in dictating the short- 
range structure of dense atomic liquids 44 . In drawing this 
analogy perhaps we should not be too surprised by the 
current observations. The reduced densities of spheres in 
the CSSM and CLJM systems are equal to 7p* ~ 0.7 
at which packing effects and short-range correlations are 
particularly pronounced (recall that the triple-point den- 
sity for the Lennard- Jones system is in the region of 0.85). 



IV. CONCLUSIONS 

In this work the structure, phase behavior, and flexo- 
electricity of model bent-core molecules have been stud- 
ied using MC computer simulations. The molecular bent 
core consists of a 'V shaped rigid array of soft spheres, 
with a transverse point dipole moment aligned along the 
C2 symmetry axis. 

In a linear seven-sphere non-polar model, isotropic, 
nematic, smectic-A, and tilted smectic-B (herringbone) 
phases are observed. With an opening angle of 160° the 
smectic A is absent, while an opening angle of 140° gives 
rise to a direct tilted smectic B - isotropic transition. The 
effects of dipolar interactions were seen to depend on the 
opening angle. In the linear-molecule systems these inter- 
actions appeared to destabilize the nematic and smectic- 
A phases. In the bent-core systems, dipolar interactions 
reduced the degree of molecular tilt in the tilted smec- 
tic phases; this is a significant observation, particularly 
since dipolar interactions have often been cited as the 
cause of spontaneous chiral symmetry breaking in some 
bent-core liquid crystals^. Chirality arises in this case 
from the correlation between molecular tilt within smec- 
tic layers and the long-range ordering of smectic-layer po- 
larizations. Although spontaneous (anti)ferroelectric or- 
der will be favored by long-range dipolar interactions, the 



9 



TABLE II: Results from NVT simulations of CSSM and CLJM systems in the nematic phase at reduced density p* and 
reduced temperature T*. 7 is the molecular bend angle, S is the nematic order parameter, el — eia 2 /p and e% — a 2 / p, are 
the reduced splay and bend flexoelectric coefficients, respectively, and P a is the average of the a component of the system 
polarization. Digits in brackets denote the estimated statistical uncertainty in the last figure based on two standard deviations. 



7 / degrees 


S 


e* 


es 




Py/H 


P./H 






CSSM, T* 


= 3, p* -- 


= 0.0888 









0.877 


0.13(4) 


0.07(3) 


-0.21(4) 


0.77(6) 


1.43(2) 


5 


0.861 


-0.01(4) 


0.13(4) 


-0.31(6) 


-0.38(6) 


0.02(2) 


10 


0.844 


0.16(4) 


0.03(4) 


0.27(6) 


-0.38(6) 


-0.10(2) 


15 


0.807 


-0.00(4) 


0.06(5) 


0.90(6) 


-1.04(6) 


0.56(2) 


20 


0.749 


0.37(5) 


-0.17(6) 


0.09(6) 


—0.55(6) 


0.67(4) 


25 


0.601 


0.51(7) 


-0.31(7) 


-1.10(6) 


1.00(6) 


-0.75(4) 






CLJM, T* 


= 6.5, p* 


= 0.1155 









0.954 


0.2(1) 


0.04(3) 


-0.13(6) 


-0.78(6) 


-0.13(4) 


5 


0.943 


-0.07(4) 


0.04(4) 


0.09(6) 


-0.42(6) 


0.13(2) 


10 


0.909 


-0.07(5) 


0.09(6) 


-0.29(6) 


-0.17(6) 


0.06(2) 


15 


0.867 


0.02(6) 


-0.10(7) 


-0.18(6) 


-1.06(6) 


0.10(2) 


20 


0.813 


-0.50(7) 


0.01(8) 


-2.03(6) 


-0.04(6) 


0.16(2) 


25 


0.739 


0.50(9) 


0.6(1) 


2.41(6) 


0.53(6) 


-0.72(2) 


30 


0.466 


0.40(1) 


0.5(1) 


-1.04(6) 


2.93(6) 


0.56(4) 



observation that these same interactions reduce molecu- 
lar tilt seems to suggest that there may be another ex- 
planation for chirality in banana liquid crystals^l. 

Real bent-core liquid crystals often possess flexible tail 
groups at the ends of the rigid bent core, and so CBMC 
simulations of a flexible-rigid-fiexible model were per- 
formed. Each molecule consisted of a five-sphere rigid 
bent core, and a three-segment flexible tail attached to 
each end. We could only find smectic and isotropic 
phases, but significantly the smectic phases showed no 
spontaneous tilt. This is probably due to the tails pro- 
viding a lubricating barrier between the smectic layers 
that serves to decorrelate the order within neighboring 
layers. Hence, the addition of molecular flexibility is 
likely to mitigate against the type of entropic 'sawtooth' 
mechanisms that have been found to stabilize antiferro- 
electric ordering in hard-particle bent-core models 6 . The 
addition of the flexible tails was also seen to give rise to 
significant spatial fluctuations in the smectic layers. It is 
worth pointing out that the smectic phases of rigid model 
molecules are often far more ordered than real smectics 
(as evidenced by scattering experiments^). The intro- 
duction of molecular flexibility therefore brings the model 
systems in to better correspondence with experiment. 

Finally, the flexoelectric properties of non-polar seven- 
sphere bent-core molecules - with and without attractive 
interactions - have been studied by calculating the splay 
and bend coefficients in the nematic phase using fluctu- 
ation relations derived from linear-response theory^S^. 
An immense investment of computational effort was re- 
quired to obtain reasonable results via this route, which 



serves to highlight how careful one must be in evaluat- 
ing the required formulae. Nonetheless, our results show 
that a significant flexoelectric response can be measured 
for opening angles below about 150°. With typical molec- 
ular dimensions and dipole moments, the measured flexo- 
electric coefficients are in the region of 10 pC m _1 which 
is in excellent agreement with experiment. The flexoelec- 
tric response of real bent-core liquid crystals is often at- 
tributed largely to dipolar and quadrupolar interactions, 
but our results show that the molecular shape is also sig- 
nificant. This shouldn't be too much of a surprise, since 
most thermotropic nematics are, after all, dense molec- 
ular liquids, and it is well known that the structure and 
dynamics in such systems are dictated by short-range re- 
pulsive interactions. We are in no way suggesting that 
electrostatic interactions are insignificant, and we have 
not studied fiexoelectricity in dipolar or quadrupolar sys- 
tems because of the computational effort which will prob- 
ably be required to obtain reliable results. A systematic 
study of this point is required, and will hopefully be the 
subject of future papers. 
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